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Abstract 

We discuss the flavor conversion of supernova neutrinos in the three-flavor 
mixing scheme of neutrinos. We point out that by neutrino observation from 
supernova one can discriminate the inverted hierarchy of neutrino masses from 
the normal one if s 2 3 > a few x 10 , irrespective of which oscillation solution 
to the solar neutrino problem is realized in nature. We perform an analysis 
of data of SN1987A and obtain a strong indication that the inverted mass 
hierarchy is disfavored unless s| 3 < a few x 10 -4 . 
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The phenomenon of neutrino flavor transformation which was discovered in atmospheric 
neutrino observation by Kamiokande and Superkamiokande experiments |l| implies new 
neutrino properties, tiny masses and unexpectedly large mixing angles. It constitutes, at 
present, the unique evidence for physics beyond the standard model of particle physics. The 
persistent discrepancy between the observed and the calculated flux of solar neutrinos 
provide another evidence in favor of these new neutrino properties. It is worth to note that 
these robust evidences perfectly fit into the standard three-flavor mixing scheme of neutrinos. 
Recently it is reported that the third member, the r neutrinos, are experimentally detected 
for the first time ||. 

Yet, there still remain many unknowns in the properties of neutrinos. Among other 
things we have no experimental clue on the question of neutrino mass pattern, i.e., the 
normal vs. the inverted mass hierarchies, either from any terrestrial experiments or solar and 
atmospheric neutrino observation. Here, we mean, by normal and inverted mass hierarchies, 
the mass pattern m 3 ^> mi ~ m 2 and mi ~ m 2 ^> m 3 , respectively. In our notation 
Am^j = m? — mf, and Am\ 2 and Amf 3 ~ Am^ are the mass squared differences that are 
related with the solar and the atmospheric neutrino oscillations, respectively. 

The issue of normal vs. inverted mass hierarchies is related to the deep question of under 
what discipline nature organizes the neutrino mass spectrum. If the inverted hierarchy is 
the case the discipline is clearly quite different from that used to organize the quark sector. 
Determination of which mass hierarchy is realized should imply powerful constraints on 
model building of neutrino masses and mixing. 

Experimentally the problem of normal vs. inverted mass hierarchies is one of the key 
issue regarding the question of whether one can measure the absolute neutrino masses. The 
only imaginable way of detecting neutrino masses of the order of J Aml tm ~ (0.04 — 0.07) 
eV is the neutrinoless double beta decay experiments, which is of course possible only for 
Majorana neutrinos. If the normal mass hierarchy is the case, one would have to go down to 
an order of magnitude smaller value of < m v >, the observable parameter in the experiments, 
to ~ O.OOleV |§§. 

In this paper we point out that neutrino observation from a galactic supernova can judge 
whether the normal or the inverted mass hierarchies is realized in nature, unless the mixing 
angle #13 is extremely small, sf 3 < a few x 10~ 4 . Moreover we perform an analysis of the 
data of SN1987A in Kamiokande and 1MB detectors || and obtain a strong indication that 
it disfavors the inverted mass hierarchy. 

Of course, this conclusion must be checked against the direct determination of the sign 
of Amf 3 which will be done in future long-baseline accelerator experiments [jJHH- However, 
the result we describe in this paper appears to be the unique hint which is available before 
such experiments are actually done. 



2 



We start by summarizing the common knowledges on neutrinos from supernova (SN) 
|T0| and their properties inside neutrinosphere [p]-|l3|. 

(1) Consideration of energetics of SN collapse indicates that almost all (~ 99%) of the 
gravitational binding energy of neutron star is radiated away via neutrino emission. The 
total energy is estimated to be several xlO 53 erg, and it is expected that the energy is 
equipartitioned into three flavors in a good approximation fil]fl4]j . 

(2) It is discussed that the shape of the energy spectra of various flavor neutrinos can be de- 



scribed by a "pinched" Fermi-Dirac distribution ||15|| . The pinched form may be parametrized 
by introducing an effective "chemical potential". 

(3) There is no physical distinction between and v T and their antiparticles in neutri- 
nosphere. It is because and are not energetic enough to produce muons by the charged 
current interactions, and the neutral current cross sections of v and v are similar in magni- 
tude. Therefore, we collectively denote them as "heavy neutrinos" in this paper 

(4) The location of neutrinosphere of heavy neutrinos, and u T , is believed to be in deeper 
place than v e and v e in SN. It is due to the fact that the heavy neutrinos have weaker 
interactions with surrounding matter; they interact with matter only via the weak neutral 
current, whereas, v e and v e do have additional charged current interactions. Hence, the 
heavy neutrinos have to have deeper neutrinosphere because their trapping requires higher 
matter density compared to those required for v e and v e . 

This last feature is of crucial importance for our business. It implies that the heavy 
neutrinos are more energetic when they are radiated off at neutrinosphere because the tem- 
perature is higher in denser region. It may be characterized by the temperature ratios of v e 
and v e to v heavy 

T = ~ ~ 1.4-2.0 (1) 

J-U e -L&e 

according to the simulation of supernova dynamics which is carried out in Ref. [TT|-|T3[. 
Despite possible slight difference it should be a reasonable approximation of ignoring tem- 
perature difference between u h and v h . 

We now turn to the the neutrino flavor conversion in supernova (SN), the core matter in 
our discussion in this paper. In fact, it has a number of characteristic features which makes 
SN unique among other astrophysical and terrestrial sources. 

(i) Because of extremely high matter density inside neutrinosphere all the neutrinos with 



cosmologically interesting mass range, m u < 100 eV, are affected by the MSW effect |T6 



*The terminology implicitly assumes that the normal mass hierarchy is the case. Nevertheless, 
we will use it even when we discuss the inverted mass hierarchy. 



3 



(Earlier references on the MSW effect in supernova include Ref. ||17|| .) Consequently, the 
three neutrino and three antineutrino eigenstates have two level crossings, first at higher (H) 
density and the second at lower (L) density, inside SN as schematically indicated in Fig. 1. 

(ii) The key question in the neutrino flavor conversion in SN is whether the H level crossing 
is adiabatic or not. If it is adiabatic, then the physical properties of neutrino conversion is 
simply u e — Vheavy exchange in the normal mass hierarchy. It should be emphasized that this 
feature holds irrespective of the possible complexity of the solar neutrino conversion which 
governs the L resonance. These key features have been pointed out in our earlier paper, Ref. 

(iii) The second important question is if the neutrino mass spectrum adopts the normal or 
inverted mass hierarchies. If the mass hierarchies is of normal (inverted) type, the H level 
crossing is in the neutrino (antineutrino) channel. 

For a recent comprehensive treatment of neutrino flavor conversion in SN in the frame- 
work of three-flavor mixing, see Ref. [|T^]. 

The last two remarks are crucial in our business. It will allow us to determine which mass 
hierarchy is realized by analyzing neutrino data from SN without knowing the parameters 
in the solar neutrino solution. Notice that this statement is valid not only for the MSW but 
also for the vacuum solar neutrino solutions. 

Before going on, let us elaborate (ii) because it is of crucial importance in our argument. 
The three solid lines at the right most end of the level crossing diagram, Fig. la, represent 
three neutrino mass eigenstates in matter at the core of SN. The neutrino evolution in SN 
obeys the Schredinger-like equation 
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where a[x) = \^2GfN £ (x) indicates the index of refraction with Gf and N e (x) being the 
Fermi constant and the electron number density, respectively. In (fj) U denotes the leptonic 
flavor mixing matrix, the Maki-Nakagawa-Sakata (MNS) matrix pO) . 

Because of an extreme density near the neutrinosphere, > 10 10 g/cm 3 , the matter term 
diag[a(a;), 0, 0] = Hq must be dominant over the other term. One can then formulate 
perturbation theory in which one takes the matter term H Q as unperturbed part and the 
other one as perturbation. It is a degenerate perturbation theory and one has to diagonalize 
the 2x2 subspace to obtain the zeroth-order wave functions and the first-order correction 
to the energy eigenvalues. The resulting three matter-mass eigenstates at neutrinosphere in 
the case of normal mass hierarchy are as follows: 
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z/g 1 is almost pure u e . v™ and v™ are certain superposition of heavy neutrinos, and u T , 
but with negligible v e component. 

The left end of Fig. la describes the three vacuum mass eigenstates V\, and u 3 . 
Therefore, if the H resonance is adiabatic, the z/3 state has the same physical properties, 
e.g., temperature, as v e inside neutrinosphere. By unitarity, the properties of v\ and 
states in vacuum, that is when they get out from the star, must be determined equally 
accurately by physical properties of and v T inside neutrinosphere. Simce there is no 
physical distinction between and v T there, the properties of v\ and are not affected by 
the nature of the L level crossing, adiabatic, moderately nonadiabatic, or nonadiabatic. We 
suspect that this feature remains true even when the two level crossings come close so that 
the independent two-resonance approximation may not be completely valid. 

This completes the argument to show that if the H resonance is adiabatic the net effect 
of the neutrino flavor conversion in SN is simply the v e - v heavy exchange, irrespective of the 
nature of the L level crossing.^] This is nothing but the feature that is called as n v e -v T 
exchange" in Ref. [18| . 



The adiabaticity of the H resonance is guaranteed if the following adiabaticity parameter 
7 is much larger than unity at the resonance point: 



7= 



Am 2 sin 2 26 
2E cos 26 
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Here, we assumed that the density profile of the relevant region of the star can be described 
as p(r) = Po( r / r o) -n to obtain the second line in the above equation, where r = 6.96 x 10 10 



cm denotes the solar radius. With the choice n = 3 and po ~ 0.1 g/cc [21 

2/3 



, we get, 



7 ~ 0.63 x 



sin 2 6 



13 



10- 



Am 



10" 3 eV 2 



E 



-1 -2/3 



20 MeV 



(4) 



for the small value of #13. Since the conversion probability P is approximately given by 
P = exp[— 17], s 2 3 > a few x 10 -4 assures adiabaticity in a good accuracy. 

Now we notice that the basic elements for the argument toward disfavoring inverted 
mass hierarchy is actually very simple. Because of (iii), the resonance is in the antineutrino 
channel if the inverted mass hierarchy is the illustrated in Fig. lb. It means that, 



^One can of course confirm this result by putting Ph = in the relevant equations in Ref. [19], 
where Ph denotes the nonadiabatic "jump" probability. 
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if the H resonance is adiabatic, all the z/ e 's at neutrinosphere are converted into heavy 
antineutrino states, and vice versa. It is also known that if H resonance is adiabatic, final 
v e spectrum at the detector is not affected by the earth matter effect jl9fl.fl 

Since the z/ e -induced charged current reaction is dominant in water Cherenkov detector, 
one can severely constrain the scenario of inverted mass hierarchy by utilizing this feature 
of neutrino flavor transformation in SN. When the next supernova event comes it can be 
used to make clear judgement on whether the inverted mass hierarchy is realized in nature, 
a completely independent information from those that will be obtained by the long-baseline 
neutrino oscillation experiments. 

While waiting for the next galactic SN, let us perform an analysis of the data of neutrinos 
from SN1987A to gain a hint to the problem of the mass pattern which we want to solve. In 
the following analysis, we assume that Si 3 is large enough, sf 3 > a few x 1CT 4 , to guarantee 
the adiabaticity of the H resonance. 

In fact, very similar analyses have been done by several authors P3| , |2"3]| . We may be able 
to characterize our work in comparison with theirs in the following way; We formulate the 
problem in a proper setting of the three-flavor mixing scheme of neutrinos, which is essential 
for the SN neutrinos. In due course, we will try to clarify how conclusions obtained in earlier 
works are to be interpreted, or to be conditioned in the three-flavor framework. 



We follow Jegerlehner, Neubig and Raffelt |24j] who employed the method of maximum 
likelihood. We define the Likelihood function as follows |2"| : 



/ POO \ -^obs 

C = C exp ( - / n(E)dE J[ n(E t ), (5) 



i=l 



where iV bs is the total number of experimentally observed events and the C is some con- 
stant which is irrelevant for our purpose of parameter estimation and the determination of 
confidence regions. Here, n(E) is the expected positron energy spectrum at Kamiokande or 
1MB detector which is computed taking into account the detector efficiency as well as energy 



resolution in the same way as in Ref. [24]. For a combined analysis of the Kamiokande and 



^The reason is as follows. Let us first note that state which carry the original v e spectrum 
oscillate very little into v e in the earth because \Am1 3 \/E is much larger than the earth matter 
potential and also because #13 is small p2[ . Lherefore, the oscillation in the earth takes place 
essentially only between v\ and U21 decoupling the state. It would lead to regeneration of v e but 
it would not give any significant effect for the v e component at the detector because both v\ and 
z?2 carry original energy spectrum of heavy flavors at the neutrinosphere. 
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1MB detectors, the likelihood function is defined as the product of the likelihood function 
for each detector. 

We draw in Fig. 2 equal likelihood contours as a function of the heavy to light tem- 
perature ratio t on the space spanned by v e temperature and total neutrino luminosity by 
giving the neutrino events from SN1987A observed by Kamiokande and 1MB detectors 
In addition to it we introduce an extra parameter r\ defined by L Vx = L Px = rjL Ve = rjL Pe 
which describe the departure from equipartition of energies to three neutrino species and 
examine the sensitivity of our conclusion against the change in the SN neutrino spectrum. 



For simplicity, as in Ref. |24|, we set the "effective" chemical potential equal to zero in the 
neutrino distribution functions because we believe that our results would not depend much 
even if we introduce some non-zero chemical potential. 

At t = 1, that is at equal v e and v e temperatures, the 95 % likelihood contour marginally 
overlaps with the theoretical expectation [13[] represented by the shadowed box in Fig. 2. 



When the temperature ratio r is varied from unity to 2 the likelihood contour moves to the 
left, indicating less and less consistency, as r increases, between the standard theoretical 
expectation and the observed feature of the neutrino events after the MSW effect in SN is 
taken into account. This is simply because the observed energy spectrum of v e must be 
interpreted as that of the original one of Vheavy, in the presence of the MSW effect in the 
anti-neutrino chanel, which implies that the original v e temperature must be lower by a 
factor t than the observed one, leading to stronger inconsistency at larger r. 

The solid lines in Fig. 2 are for the case of equipartition of energy into three flavors, 
r) — 1, whereas the dotted and the dashed lines are for r\ = 0.7 and 1.3, respectively. We 
observe that our result is very insensitive against the change in rj. 

We conclude that if the temperature ratio r is in the range 1.4-2.0 as the SN simulations 
indicate, the inverted hierarchy of neutrino masses is disfavored by the neutrino data of 
SN1987A unless the H resonance is nonadiabatic.| 

For completeness, we summarize the features of neutrino events that we expect in the 
three-flavor mixing scheme of neutrinos. They differ depending upon the normal or the 
inverted mass hierarchies and on the nature of the H and the L level crossings. 
The case of normal mass hierarchy 

There is no level crossing in the antineutrino channel, apart from the possibility that the 



solar mixing angle is in the "dark side" in the parameter space |26|, 612 > tt/4, which we 



do not consider in this paper. Therefore, as far as #13 is not large, which is indicated by 



CHOOZ result p2[, third state u 3 is essentially decoupled from the other antineutrinos and 



§This result has been announced at two international conferences, Dark2000 and NOW2000 [25]. 
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as far as v e signal is concerned, the problem can be approximately reduced to the two flavor 
mixing scheme which is well explored by the previous works p3|j2^| . 

The conclusion reached can be summarized as follows: If the small mixing angle (SMA) 
MSW is the solution to the solar neutrino problem there is only a minor effect because 
neither the vacuum oscillation nor the earth matter effects are effective because of small 612. 
If the large mixing angle (LMA) or low Am 2 (LOW) MSW or vacuum oscillation (VO) is 
the solution, we have a potential trouble because a good fraction of v e is transformed into 
Vheavy and vice versa. One can repeat the similar analysis as we did for Fig. 2 and would 
conclude that all of these solutions with large mixing angle are disfavored |23], though less 
convincingly than the case of inverted mass hierarchy. 

Fortunately, the earth matter effect helps us to cure the trouble at least to some extent 
for the case of LMA MSW solution [23], 24]. We present in Fig. 3 the result of the analysis 
using the same likelihood method. We employ a particular set of parameters of the LMA 
MSW solution and compare the behavior of the likelihood contours with and without earth 
matter effect. For simplicity, we set #13 = but our result does not change much as long as 



the parameter is under the CHOOZ bound ||22[| . In the present analysis we only deal with 
the case of equipartition of energy, 7] = 1, because we already know that the results are 
not sensitive to rj. As one observes in Fig. 3, the earth matter effect cures the discrepancy 
between the likelihood contours and the theoretical expectation to some extent, but not 
completely. 

For the case of LOW MSW as well as VO solutions, there is no significant earth matter 



effect and the results are essentially the same as was presented in Fig. 5 in Ref. ||24|| , and 
hence we show no plot for these solutions. 

There exist interesting effects in the neutrino channel because they have level crossings. 
We have to discuss three cases separately; the H resonance is (a) adiabatic, (b)moderately 
nonadiabatic, and (c) nonadiabatic. 

In the case (a) the net effect is v e -Vh eavy exchange independent of the nature of the L 
resonance, as we have discussed extensively. The characteristic signature of harder spectrum 
of v e which comes from Vheavy i n neutrinosphere is: 

(i) Enhancement of forward peaking elastic scattering events at high energies which should 
be observable in water Cherenkov detectors f!8 |, and 

(ii) Enhanced oxygen-induced events due to a steep rise of the cross section at energies higher 



than > 30 MeV p7[ , which could be separated from the dominant isotropic z7 e absorption 
events due to a moderate backward peaking of the events |f28|,f2"9|j . 

In the case (c) one can disregard the third state and the problem is the pure two-flavor 
fe-fheavy transformation. Then, if the LMA or LOW MSW is the solution to the solar 
neutrino problem, there is a significant conversion with similar experimental signatures as 
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the case (a). We expect that nothing happens in SN in the case of VO solution, and about 
equal mixture of u e and v heavy (for sin 2 2#i 2 — 1) would result due to vacuum oscillations 
which would also lead similar but somewhat weaker experimental signatures as the case (a). 
In the case of SMA MSW solution the conclusion depends upon the profile of the outer 
envelope of progenitor star and it is difficult to draw definitive conclusions. 

In the case (b) in which the conversion probability at the H resonance is, say, ~ 0.5, then 
one can expect some similar but weaker effects with the case (a) independent of the solar 
neutrino solution, and in addition the effects which depend upon the solar solutions. 
The case of inverted mass hierarchy 

The structure of the level crossing is very simple in this case as shown in Fig. lb. 
There is the H resonance in the antineutrino channel and the L resonance is in the neutrino 
channel, since we do not consider the solar parameter which is in the "dark side". Then, 
the predictions to the neutrinos are exactly the same as the ones for case (c) in the normal 
mass hierarchy. 

In the antineutrino channel there is the H resonance and we have discussed extensively 
what happens if it is adiabatic, the ). So we only need to discuss the case (b) and (c). 

If the H resonance is nonadiabatic, the case (c), the effect of neutrino mixing is the same as 
in the case of normal mass hierarchy. In the case of (b), moderately nonadiabatic case, v e - 
Vheavy transformation occurs with the probability 1 — Ph, and it would imply the similar but 
milder effect than that we have obtained with adiabaticity of the H resonance. If the next 
galactic supernova is detected by Superkamiokande, then we will be able to discriminate the 
moderately nonadiabatic case from the adiabatic one. 

After completion of our work we became aware of the paper by Lunardini and Smirnov 
|30| in which some points related with our work are mentioned, but with particular emphasis 
primarily on the detector-dependent earth matter effect. 
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FIGURES 




FIG. 1. The schematic level crossing diagram for the case of (a) normal and (b) inverted mass 
hierarchies considered in this work. The circles with the symbol H and L correspond to resonance 
which occur at higher and lower density, respectively. 
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FIG. 2. Contours of constant likelihood which correspond to 95.4 % confidence regionsfor 
the inverted mass hierarchy under the assumption of adiabatic H resonance. From left to right, 
r = T Px /T Pe = T Vx /T Ve = 2,1.8,1.6,1.4,1.2 and 1.0 where x = n,r. Best-fit points for T Pe and 
Ef, are also shown by the open circles. The parameter rj parametrizes the departure from the 
equipartition of energy, L Vx = L Px = rjL Ve = 7]L Pe (x = //, r), and the dotted lines (with best fit 
indicated by open squares) and the dashed lines (with best fit indicated by stars) are for the cases 
rj = 0.7 and 1.3, respectively. Theoretical predictions from supernova models are indicated by the 
shadowed box. 
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FIG. 3. Contours of constant likelihood corresponding to 95.4 % C.L. for the large mixing an; 
MSW solution (a) without and (b) with earth matter effect. We have taken mixing parameters 
Am 2 = 3 x 1(T 5 eV 2 and sin 2 26 = 0.8 for LMA MSW solution. 
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